Melting of a frustration- induced dimer crystal and incommensurability in the J1-J2 two-leg ladder 
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The phase diagram of an antiferromagnetic ladder with frustrating next-nearest neighbor couplings along the 
legs is determined using numerical methods (exact diagonalization and density-matrix renormalization group) 
supplemented by strong-coupling and mean-field analysis. Interestingly, this model displays remarkable fea- 
tures, bridging the physics of the J1-J2 chain and of the unfrustated ladder. The phase diagram as a function 
of the transverse coupling Jx and the frustration J2 exhibits an Ising transition between a columnar phase 
of dimers and the usual rung-singlet phase of two-leg ladders. The transition is driven by resonating valence 
bond fluctuations in the singlet sector while the triplet spin gap remains finite across the transition. In addition, 
frustration brings incommensurability in the real-space spin correlation functions, the onset of which evolves 
smoothly from the J1-J2 chain value to zero in the large- Jx limit. The onset of incommensurability in the 
spin structure-factor and in the dispersion relation is also analyzed. The physics of the frustrated rung-singlet 
phase is well understood using perturbative expansions and mean-field theories in the large- Jx limit. Lastly, we 
discuss the effect of the non-trivial magnon dispersion relation on the thermodynamic al properties of the sys- 
tem. The relation of this model and its physics to experimental observations on compounds which are currently 
investigated, such as BiCu2P06, is eventually addressed. 

PACS numbers: 75.10.Kt, 75.40.Mg, 75.10.Jm, 75.10.Pq 



Ladder materials offer a unique playground to improve our 
understanding of the subtleties arising from quantum effects 
in low dimensional geometries.' For antiferromagnetic 5* = 
1/2 Heisenberg models, it is well-known that ladders with an 
even number of legs ni display short-range correlations for 
any non-zero inter-chain coupling J±^ and a finite spin gap 
As ~ J^ cxp(— ajif), whereas the gapless quasi-long-range 
ordered state of a single S = 1/2 chain is robust for ladders 
having an odd number of legs. ^ Ladder physics has been inten- 
sively explored during the last two decades both theoretically 
and experimentally, in particular regarding spin gap physics, 
impurity effects,^ field-induced magnetization processes,"* su- 
perconductivity in hole doped systems,' etc. Despite their 
rather simple geometry and the huge amount of studies, ladder 
systems remain a topic of current interest. Newly synthesized 
two-leg ladder materials which exhibit sizable spin gaps ~ 1 
meV have recently emerged,"*-^ thus opening the possibility to 
close the gap with an external magnetic field. 

The ground-state of a spin-1/2 Heisenberg two-leg lad- 
der is a genuine quantum state which has no classical ana- 
log. Schematically, the short-range nature of the spin correla- 
tions can be encoded in the so-called resonating valence bond 
(RVB) picture with short range pairwise singlet bonds fluctu- 
ating over a few lattice sites. ^ Such a state, sometimes called 
"rung-singlet" (RS) because the strongest antiferromagnetic 
coiTelations are along the vertical rung bonds, is character- 
ized by confined triplet excitations and a magnetic response 
to a local perturbation such as a non-magnetic impurity which 
remains confined in its vicinity, as probed by NMR. Another 
kind of spin-liquid which has also been intensively studied by 
many authors, back to Majumdar and Ghosh (MG) 40 years 
ago,^ occurs the so-called Ji — J2 chain. Used as a prototype 
for modeling the spin-Peierls transition in frustrated chains 
such as CuGeOa, it displays a Kosterlitz-Thouless transition 
at a finite frustration ratio J2/J1 from a quasi-Neel-ordered 



state to a dimerized gapped state, which spontaneously breaks 
the (discrete) translational symmetry and is two-fold degen- 
erate. The model under study in this manuscript bridges the 
physics of these two famous spin-liquids by considering two 
coupled J1-J2 chains, or equivalently, a two-leg ladder with 
frustrating couplings along the legs (see Fig. 1). It realizes a 
transition from a columnar dimer phase to the RS phase. 

The search for a transition between the RS phase and dimer- 
ized phases in spin-1/2 two-leg ladders was initially triggered 
by Ref. 8 which proposed two scenarios. Based on the rewrit- 
ing of the low-energy physics in terms of four gapped Ma- 
jorana fermions fields^ comprising one in the singlet sector 
and three in the triplet sector Taking into account a four 
spin interaction** leads to closing the triplet gap or the singlet 
gap while the other sector remains gapped. The first scenario 
was realized in the well-studied two-leg ladder model with 
ring exchange,'" displaying a transition to staggered dimer 
phase. The second scenario has been realized following the 
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FIG. 1. The J1-J2 two-leg ladder model (top). Qualitative sketches 
of the two phases realized : the columnar dimer phase (middle) and 
the rung-singlet phase (down). 



initial proposal by taking into account four spin interaction 
with a negative coupling." Restricting models to two-body 
exchanges only, the main frustrated models which have been 
investigated to stabilize dimerized phases are two-leg ladders 
with cross-coupled, zig-zag and inequivalent exchanges on 
plaquettes.'^ Other models considering explicit dimerization 
can naturally yield staggered and columnar dimer orders.'^ 
In the context of the confinement-deconfinement transition,'^ 
the J1-J2 ladder model for which there is an experimental re- 
alization,'* was studied using the low-energy physics'^ and 
numerics, yet limited to the MG point'-'"^ J2I J\ = 0.5 and 
with little evidence. 

Here, we present an exhaustive study of the full phase dia- 
gram of the Ji- J2 laddermodel of Fig. 1, motivatedby several 
aspects such as the shape of the whole phase diagram, the na- 
ture of the excitations, and the evolution of the incommensu- 
rability upon increasing frustration. Besides the mere theoret- 
ical interest, this model is also relevant to the spin ladder ma- 
terial BiCu2P06 where several studies have pointed towards a 
scenario involving a sizable intra-chain frustration J2/ Ji-^'"^^ 
In particular, inelastic neutron scattering experiments display 
incommensurate spin correlations which originate from mag- 
netic frustration.^" Based on numerical and analytical calcu- 
lations, we report a careful study of the incommensurate re- 
sponse upon increasing the frustrating coupling. 

The rest of the paper is organized as follows. We first 
present in section I the full phase diagram of the model which 
display two gapped phases (columnar dimer and rung-singlet) 
together with a crossover line which signals the onset of in- 
commensurability. Numerical results are presented in sec- 
tion II where the nature of the quantum phase transition, low 
energy spectra and incommensurate correlations are studied. 
In section III, we present analytical results based on strong- 
coupling expansions at large Jx and bond-operator mean-field 
approximation for the triplet excitation branch. They account 
for the evolution of the spin gap and incommensurate wave- 
vectors in the RS phase obtained by simulations. The temper- 
ature dependence of the magnetic susceptibility is presented 
in section IV. Finally, section V A concludes the paper and 
discusses some experimental issues and open questions. 



I. MODEL AND PHASE DIAGRAM 

We study the following Hamiltonian on a spin- 1/2 two-leg 

ladder and depicted on figure 1: H = X]j=i ^i^ ^^^ the three 
couplings 

hi — J\ [Si^i • Si-|_i,i + Si^2 • Si+1,2] (1) 

+J2 [Si,i • Si+2,1 + Si_2 ■ Si+2,2] (2) 

+Ji So • S,,2 , (3) 

where S^j is the spin-1/2 operator acting at site i of leg 
j and the Js are the magnitude of the various couplings 
which are here taken to be antiferromagnetic (J > 0). The 
zero-temperature phase diagram is derived using numerical 
methods well suited for frustrated systems : exact (or Lanc- 
zos) diagonalization (ED) and density-matrix renormalization 



group^^ (DMRG), keeping up to 2000 states in each block. 
Calculations are carried out using periodic boundary condi- 
tions (PBC) for ED and open boundary conditions (OBC) for 
DMRG. We now describe below the main features of the ob- 
tained phase diagram. 

This model includes two limiting cases of two different 
spin-liquid phases which have been extensively studied. We 
summarize below the main known results about these case: 

(i) When Jj^ = 0, the model corresponds to two decoupled 
J\-J-2. chains. A single J1-J2 chain displays a transition 
between the gapless phase of the Heisenberg chain to a 
gapful dimerized phase for-^ J2,c — 0.24117 Ji. For 
the special Majumdar-Ghosh (MG) point^ J2 = 0.5 Ji, 
the ground-state can be written in a simple product of 
dimers. The gap and the spin correlation length can be 
computed exactly at this point. Beyond this remarkable 
mathematical feature, the MG point physically turns 
out to be the point at which the correlation length is 
minimal and beyond which the real-space spin correla- 
tions start to be incommensurate. As incommensurabil- 
ity will be discussed in details in this paper, we already 
give the three definitions and notations of wave-vectors 
signaling incommensurability in the three different ob- 
servables: 



real-space spin correlations 
maximum of the spin structure-factor 
minimum of the magnon dispersion relation 



(4) 



In the J1-J2 chain, due to the finite correlation length, 
the spin structure-factor (Fourier transform of the real- 
space correlations) has a maximum at an incommensu- 
rate wave-vector q above the so-called Lifshitz point^"* 
Ji.L — 0.5206 J\. Notice that the onset of incommen- 
surability in the minimum of the triplet dispersion re- 
lation at q* is expected to happen at a third different 
value of J2/J1 (as for the analogous spin-1 chain-^), 
but remains hard to study numerically because of lim- 
ited sizes with exact diagonalization. In the large J2I J\ 
limit of two weakly-coupled chains, the behaviors of the 
spin gap and q follow reasonably well the bosonization 
predictions.-''-^^ 

(ii) When J2 = 0, the model is that of two coupled Heisen- 
berg chains, i.e. the Heisenberg spin ladder. A spin gap 
in the magnetic excitation opens linearly with the trans- 



verse coupling Jx with logarithmic corrections • at 
small J^x, leading to the RS phase. The naming orig- 
inates from the large- J^ regime in which the picture 
of the ground-state boils down to dimers on each rung. 
However, the correct physical picture of the isotropic 
ladder is that of an RVB spin-liquid* rather than that of 
pinned dimers along rungs. 

The main result of the paper is the phase diagram of Fig. 2 
which displays two main phases (apart from the critical line 
,/j^ = and J2 < J2,c) : a columnar phase of dimer 
which breaks the translational symmetry and the RS phase. 
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FIG. 2. (color online). Phase diagram of the J1-J2 ladder. The 
commensurate-incommensurate regimes separated by the dashed 
green line corresponds to the incommensurability emerging in the 
real-space spin correlation functions. The violet transition line at fi- 
nite Jx is of the Ising type. The dashed violet line is the conjectured 
behavior of the transition line close to J2,c ~ 0.241 Ji. 



The columnar phase, first proposed by Vekua and Honecker'^ 
along the MG line J2 = 0.5 Ji, emerges naturally as the 
combination of two dimerized chains with their dimer pat- 
terns in phase. The physical picture of the transition is the 
following : starting from J± = 0, we have two dimerized 
chains between which the interchain antiferromagnetic cou- 
pling will favor RVB fluctuations leading to the alignment 
of their dimers (reminiscent of two-dimensional'^ colum- 
nar phases) and a lowering of the order. RVB fluctuations 
enhanced by larger J± eventually destroy the dimer order 
through an Ising second-order transition to the RS phase (as 
expected for a single-component order parameter, the transi- 
tion coiTesponding to breaking a discrete Z2 symmetry). As 
we will see, the triplet gap remains finite across the phase di- 
agram and on the transition line between these two distinct 
spin-liquid phases. This melting of the dimer crystal is driven 
by the low-energy fluctuations which lie in the singlet sector 

The transition line has a remarkable non-monotonic behav- 
ior showing a reentrance of the RS phase at small J± when 
increasing the frustration J2. This shape can be understood 
on simple qualitative grounds : the order parameter and spin 
gap in the J1-J2 chain (along the J± = line) typically have 
a steep increase above the transition point and then a slow 
decrease for large J2 . Considering J± as a perturbation, the 
larger the initial order, the larger J± is required to destroy 
it. This argument is all the more valid as the magnitude of 
J± along the transition does not exceed 0.4 Ji. When the 
order parameter is in the Ising scaling regime, we may use 
Eq. 7 as a rough estimate of the transition line which yields 
: J±^c oc [£'chmn(-^2)]*, with Dchain(-^2) the dimer Order pa- 
rameter of a single Ji — J2 chain. Furthermore, the numer- 
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FIG. 3. (color online), (a) Behavior of the order parameter across 
the transition at fixed J2 = 0.6. Extrapolations to the thermody- 
namic limit are performed using a polynomial fit in the ordered phase 
and Eq. (6) close to the transition point, (b) Scaling of the order pa- 
rameter around and at the critical showing the exponent of the Ising 
universality class. 



ics are particularly difficult when J \^ and the gaps are small, 
which corresponds to the behavior of the transition line close 
to J2.C — 0.24Ji and at large J2. There, we may conjec- 
ture that the transition line qualitatively follows the behavior 
of the J^ = order parameter and should therefore exhibit 
an exponential-like opening Jj^c « e~^'^i/*^'^^^''^'=) close to 



J2.C and an exponential tail J_l_c oc e 



-A'Ja/Ji 



at large J2I J\ 



The last remarkable feature of this phase diagram is the 
onset of incommensurability which is a typical signature of 
frustration, akeady present at the classical level. As already 
recalled in the case of the frustrated spin chain, the onset of 
incommensurability depends on choosing the correlations in 
real or reciprocal spaces, or looking at the dispersion relation. 
In Fig. 2, we chose the standard definition for the onset of 
incommensurability as the appearance of an incommensurate 
wave-vector in the real-space spin correlations. 

In the following sections, we explain how these results are 
obtained and describe in details the physics throughout the 
phase diagram using several observables. 



II. NATURE OF THE TRANSITION 

A. Order parameter 

In this section, we give numerical evidence for the Ising 
nature of the transition between the columnar and RS phases. 
The simplest way to distinguish the two phases is to compute 
the dimer order parameter defined around site i on leg j by 



Uij — Wi.j ' ^J 



+ i.i. 



- (S, 



-ij 



'1,21 



(5) 



Using DMRG with OBC, one can directly access Dij as a lo- 
cal order parameter We check that the dimers are aligned on 
each legs and not staggered, which would be the other possi- 
ble pattern on a ladder but energetically less favorable in this 



model. Therefore, we drop the leg index j in the following. 
Due to OBC, Friedel oscillations develop in the dimer order 
from the edges which decay exponentially in both the colum- 
nar and RS phases (away from the transition) and contain 
incommensurate oscillations in the incommensurate regime. 
Ignoring these oscillations to keep only the envelope of the 
Friedel oscillations, their typical behavior near the edges is 
expected to be^" 



D^ ~ D^ 



A- 



'X/i,; 



(6) 



with X the distance from the edge, Doo the value of the or- 
der parameter in the thermodynamic limit (in the bulk of the 
ladder), A a constant and ^dimer the correlation length associ- 
ated with the dimer fluctuations. The exponent a accounts for 
power-law corrections which are particularly relevant close to 
the transition. We expect Doo to be zero in the RS phase and 
finite in the columnar dimer phase. We keep the discussion 
of Cdimer for Scc. Ill B but wc Can already point out that it is 
a different length-scale from the usual spin correlation length 
Cspin obtained from spin correlations. 

The typical evolution of the order parameter £'l/2 at the 
middle of the ladder with increasing J± and fixed J2 =0.6 Ji 
is shown on Fig. 3(a). The vanishing of the order parameter 
at the transition is in qualitative agreement with the Ising pre- 
diction 



Doo{J±) OC (J_l,c - J±) 



1/8 



(7) 



in the vicinity of the critical point. In order to have a much 
more accurate determination of the critical point, as well as a 
test of the universality class, we use the finite-size scaling of 
the order parameter I?l/2 given by Eq. (6). In the quantum 
Ising universality class, the correlations of the order parame- 
ter have a decay exponent 1/4 at the critical point (for which 
Cdimei- = oo), which givcs an exponent a — 1/8 for the associ- 
ated Friedel oscillations. Fig. 3(b) shows the scaling behavior 
across the transition and the good agreement between the ex- 
ponent found numerically and the Ising one. 



B. Entanglement entropy 

Another supporting argument for the Ising nature of the 
transition is the value of the central charge c = 1/2 at the crit- 
ical point. In order to numerically extract c from the DMRG 
data, we use the scaling of the entanglement entropy S{x) 
with the half-block size x for OBC: 



S{x) 



6 



■In 



L + 1 



L + 1 



(8) 

where A, B and c are fitting parameters. This formula is based 
on a universal prediction for the entanglement entropy^' and 
subleading corrections-'^ {Sx+i i^x 1) which are here nothing 
but the local dimer order calculated numerically. In order to 
obtain a good scaling, one must precisely locate the transition 
point (using the previous approach) and use rather large sys- 
tem sizes (here up to i = 400). On Fig. 4, we observe a scal- 
ing compatible with c = 1/2, providing another independent 




FIG. 4. (color online), (a) Entanglement entropy vs. block size i 
fitted by Eq (8) for the critical point of Fig. 3. (b) Scaling of the 
entanglement entropy of a half-system with size L: The extracted 
central charge c ~ 0.46 is in good agreement with the Ising value 

c= 1/2. 



evidence of the Ising nature of the transition. We stress that, 
since the spin correlations remain short-range at the critical 
point, the logarithmic increase of the entanglement entropy is 
only driven by the critical dimer fluctuations. As the entan- 
glement entropy oscillations are related to the dimer order, its 
block-size dependence shows an incommensurate behavior as 
in Ref. 33. 



C. Excitation spectra from ED 

In this section, we discuss the nature of the low-energy ex- 
citations as a signature of the two phases. As a preamble, 
we recall that the deconfined spinons living on a dimerized 
chains must get confined as soon as the transverse coupling 
,/j^ is turned on,'^ since they must bind into a singlet state 
which size is controlled by J^. Hence, the spinons excitations 
are confined throughout the phase diagram, except along the 
Jj^ = and J 2 > J2.C line. We report in Fig. 5 the finite- 
size spectra as functions of the momentum k^ in the two par- 
ity sectors (defined by the transverse momentum ky) and for 
the singlet and triplet channels (identified using spin inversion 
symmetry), for a fixed frustration J2 = 0.6 Ji and increasing 
J±- 

Starting with a small J± (Fig. 5(a)), the spectrum has four 
nearly degenerate low-energy states. They clearly stem from 
the four possible singlet states corresponding to the combina- 
tion of two doubly degenerate chains. The degeneracy is lifted 
by J± which stabilizes the combinations associated with the 
columnar phase, while staggered configurations^"^ have a gap 
controlled by J±. Above these four low-energy singlet states, 
gapped singlet and triplet modes reminiscent of the dimerized 
chain"*^ develop around k^ = 0, it. Their near degeneracy is 
a signature of weakling confined spinons. We notice a min- 
imum in the magnon dispersion relation occurring at an in- 
commensurate wave-vector q* ~ 37r/4, using the notation of 
Eq. (4). Comparison between the two wave-vectors q and q* 
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FIG. 5. (color online), (a-d) Low-energy excitations spectra for fixed J2 = 0.6 J\ and increasing Jx. The up arrows indicate an approximate 
location of the incommensurate wave-vector q* corresponding to the minimum of the magnon relation dispersion. Full (dashed) symbols are 
for 1/ = 18 (1/ = 16). (e) scaling of relevant gaps very close to the transition point : the triplet excitations remain gapped while a gapless 
singlet mode develops around the wave-vector (tt, 0) as can be seen from the vanishing of the singlet gap at (vr — 2-k / L, 0). 



will be addressed in Sec. IIIC. 

Very close to the critical point (spectrum with J^ = 
0.37 Ji of Fig. 5(b)), we observe that the triplet excita- 
tions remain fully gapped, while a low-energy singlet mode 
emerges around (tt, 0) (the wave-vector of the nearly de- 
generate ground-state). The emergence of this mode is sig- 
naled in Fig. 5(e) by the scaling to zero of the singlet gap at 
(tt — 2tt/L, 0). On the contrary, the triplet modes clearly have 
a finite gap in the thermodynamic limit. This singlet mode is 
related to the algebraic decay of the dimer correlations, sup- 
porting the melting of the dimer crystal scenario for the tran- 
sition. This scenario seems to agree with the field theoretical 
approach and early numerical results of Ref. 13. 

Entering the RS phase by increasing further J^, the ground- 
state is now non-degenerate and the spectrum fully gapped 
with an incommensurate magnon branch as the lowest excita- 
tion. For J± = Ji (Fig. 5(c)), this branch enters the contin- 
uum of many-magnons excitations for small and large kx, but 
for large enough transverse coupling {J± ~ 3 Ji in Fig. 5(d)), 
the incommensurate magnon branch is well separated from 
the continuum (as usually found for Heisenberg ladders). This 
strong-coupling physics is well captured by the perturbative 
calculations of Sec. III. 



D. Incommensurate wave-vector from spin correlations 

The onset of incommensurability across the phase diagram 
is obtained by choosing the definition from real-space spin 
correlations, i.e. when q deviates from tt. Following the ap- 
proach of White and Affleck for the J1-J2 chain,^^ we extract 
the incommensurate wave-vector q with the following ansatz 
for the exponential decay of the spin correlator : 



((& 



V+x,l 



■'y+x,2J\^y 



)iS'i~S'2))(xcos{qx + ip) 



-a;/6p 



(9) 
where ^spin is the spin correlation length, which is also ex- 
tracted from the fit. The antisymmetric combination of the 



spin operators is chosen to get the signal corresponding to the 
magnon mode with ky = tt. The magnitude of the frustrating 
term required to drive q away from tt is reported on Fig. 2 
by the green dashed line. The line starts at the MG point 
J2/J1 = 0.5 for J± ~ and asymptotically decreases to- 
wards J2/J1 — in the large J± limit. There, an accurate de- 
scription of the onset line, as well as the comparison between 
the wave-vectors (4) is accessible (see Sec. IIIC). We lastly 
notice that for large J2/J1, q naturally reaches tt/2 since the 
model boils down to two weakly coupled ladders with a dou- 
bled unit cell. 



III. MAGNON DISPERSION RELATION 

This section is dedicated to a more detailed study of the 
magnon excitation branch and its main features : effect of 
the frustrating term, spin gap, correlations lengths and the on- 
set of incommensurability. The strong-coupling limit offers 
tractable analytical predictions that are compared with exact 
numerical results. 



A. Perturbative approaches at large J± 

1. Strong-coupling expansion 

A strong-coupling expansion can be carried out in the large- 
J± limit taking Ji and J2 as perturbations. These two terms 
lead to nearest neighbor and next-nearest neighbor hoppings 
for the magnons and the dispersion relation reads 



oj{k) = J± + Ji cos(fc) + J2 cos(2fc) + 



3J? 

4 



J^ 



Ji 



(10) 



In comparison to the standard ladder dispersion relation, frus- 
tration brings the harmonic 2fc leading to incommensurabil- 
ity. Remarkably, and contrary to the case of a single J1-J2 
chain, we here have an exact estimate for the wave-vector 



minimizing the magnon gap, which we wrote q* and which 
is, in principle, different from q. In this limit, if J2 < Ji /4, 
the minimum of the dispersion remains at g* = tt while for 
J2 > Ji/4, the minimum takes place at the incommensurate 
wave-vector 



with 



A 

4J2 



(11) 



which reaches 7r/2 in the large-J2 limit. Note that in the 
regime where Ji ^ J^, J2, the picture is that of two un- 
frustrated two-leg ladders with a doubled lattice spacing and 
couplings J2, J±, slightly coupled in a zig-zag way by the per- 
turbation Ji . The doubling of the unit cell naturally gives 
a minimum at q* = it/2 in the original lattice model. The 
threshold and incommensurate wave-vector (11) are the ones 
that actually appear as the angle of the helical structure of 
the classical version of both the J1-J2 chain and ladder mod- 
els, hence in the spin-wave calculation. While for a single 
chain and the small- J^ regime, quantum fluctuations strongly 
change the onset value and the q*{J2) function, the calcula- 
tion is here exact in the large- Jj^ regime. 



2. Bond-operator mean-field theory 

The shape of the dispersion relation can be refined using 
bond-operator mean-field (BOMF) theory following Ref. 37. 
The starting point is here again the large- Jj^ limit where 
the local Hilbert space on each rung is spanned by intro- 
ducing a singlet creation operator s\, and three triplet cre- 
ation operators t\ ^^, t\ and t\ ^. The local constraint s\si + 

X^Q^Q^^a = 1 is introduced using a Lagrange multiplier 
which identifies to a local chemical potential ji which we as- 
sume to be uniform from translation invariance. The Hamil- 
tonian terms (l)-(3) can then be rewritten using these creation 
operators as: 



i,a''^,oi\ 



Jl, 



+ Y^tlj^+l,a4+lSt + 4,a4+l,a«i+lS* + h-C-] 






^Si+2Si 



h.c. 



l^[sUi+tlaU.. 



in which i labels the rung, there is an implicit sum over a = 
x,y,z, and the dots accounts for triplet-triplet interactions that 
we neglect according to Ref. 37. The resulting Hamiltonian is 
then studied under a mean-field approximation relying on the 
fact that the RS phase is dominated by singlets on rungs in 
the large- J^ regime. We thus take s^ ~ s which yields the 
Hamiltonian (in Fourier space) governing the dynamics of the 
triplets: 



Hk 



=L 



E 



T -2 -2 I 



tk,at- 



k,a\ 



Afc = J±s'^ (Xi cos(fc) + A2 cos(2A:) j + — ^ - /i 



A. - 



-s^iXi cos(fc) + A2 cos(2A:) j 



where we introduce two small parameters Ai_2 = Ji.2/J±- 
This Hamiltonian is solved by a Bogoliubov transformation 
leading to 
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and uj{k) - ii^ - ^^fc 

and s are found by minimizing Eq. (12) which gives 
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J± s—^ Ai cos(fc) + A2 cos(2fc) 
2L ^ ^\ + di cos(fc) + d2 cos(2A:) 
r-^ 2 + di cos(fc) + 6,2 cos(2fc) 



2 4L ^ ^1 + di cos(fc) + d2 cos(2A:) ' 



with di = y7|^^ and rf2 = - 



2A2g 



The self-consistent 



relations are solved numerically to obtain the best {n,s). 
The dispersion relation, within the mean-field approximation, 
reads 



Lu{k) 



c-t 



^ ) Vl + di cos(fc) + d2 cos(2fc) . (13) 



As 6,2/ di = J2/J1, one obtains the same condition J2 = 
Ji/4 for the onset of incommensurability and the same in- 
commensurate wave-vector given by (11). Taking the expan- 
sion of the mean-field equation in the As, one recovers the 
same first-order correction as in Eq. (10) but the second order 
corrections are different : 

u}{k) = J_L [1 + Ai cos(fc) + A2 cos(2fc) + -{\\ + Xf) 



1 



(Aicos(fc) + A2Cos(2fc))^ 



In particular, the spin gap to second order is not well repro- 
duced. 



3. Comparison with exact diagonalization 

We can now compare these analytical results to exact diag- 
onalization results computed on ladders with up to L = 18 
rungs. As expected. Fig. 6 shows a good quantitative agree- 
ment with the predictions in the strong-coupling regime in 
which finite-size effects are pretty small as one can judge 
from the data. The strong-coupling prediction gives an overall 
good account of the dispersion relation, with the best accuracy 
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FIG. 6. (color online). Magnon dispersion relations for J± = 3Ji 
and J2 = O.6J1 from Lanczos diagonalization and compared to 
strong-coupling and BOMF theory. 



around tt. The BOMF result (13) needs a self-energy correc- 
tion-'^ of 0.7(Jf + J2)/ J\^ to match the correct spin gap. 
Once this is done, the BOMF provides a better description 
than the strong-coupling expansion, particularly in the vicin- 
ity of q* . Lastly, as for the well-known unfrustrated ladder, the 
picture from the strong-coupling limit remains qualitatively 
valid down to the isotropic ladder limit Ji ~ J±^ (see Fig. 5). 



B. Spin gap and correlation lengths 

1. Spin gap and incommensurability with increasing frustration at 
large and fixed J± 

A Striking feature of the dispersion relation (10) is that, 
at large- J^, the gap first increases with J2 before decreas- 
ing, hence passing through a local maximum. In the strong- 
coupling limit, the spin-gap As has the following behavior 



A, = 



r r T 3J2 + 

4 J I 



•^2 



J±-J2 



8J2 



3Jf 

4 



J^ 



J I 



for J2 < Ji/4 
for J2 > Ji/i 



(14) 

For J± = 3 Ji in Fig. 7(b), the agreement with the numerics 
is already pretty good. One may observe from Eq. (14) that 
the local maximum of the spin gap is actually realized for a 
value of J2 larger than the onset of the q* incommensurability 
wave-vector Neglecting terms in J2/ J_l, the maximum from 
(14) occurs at J2 = Ji/\/8 ~ 0.35 Ji, which agrees with 
Fig. 7(b). For the isotropic ladder with J± = Ji, the strong- 
coupling result is no more quantitative but the existence of 
a local maximum and its location above the onset of q* (not 
shown) are yet robust features (see Fig. 7(a)). We will see in 
Sec. IIIC that q* gets very different from q in the large- Jj^ 
limit so that the local maximum is actually not related to the 
true onset of incommensurability associated with q. 






0.5 



onset of q 

■ 



spin gap (DMRG) 
spin gap 
(strong-coupling) 

1/^ ■ 



,1 I ■ 

onset of q ' 



(a) J=J, 



0.5 



2 

S 1.5 

1 

0.5 




f 

I oiiHet of q 
■ ■ 

- (b) J^=3J^ 







0.5 



FIG. 7. (color online). Evolution of the triplet gap and inverse spin 
correlation length with J2 at fixed J±. The onset of the incommensu- 
rability in spin correlations (g), spin-structure factor (q) and disper- 
sion relation (q*) are given when available, (a) the isotropic ladder 
situation J± — J\, where we have q ^ q. The dashed line is just 
a guide for the eyes, (b) in the strong-coupling regime Jx = 3Ji 
where q and q are very different. The line is the strong-coupling 
behavior of Eq. (14). 



One may wonder whether the spin correlation length <^spin 
obtained from Eq. (9) has also a singular behavior as it has 
been observed for instance in the spin- 1 chain^^'^^ and whether 
its location is related to the onset of incommensurability. 
From Fig. 7, we do observe that £^^1^ displays a local singular 
maximum close to the onset of incommensurability associated 
with q (and not q* or q). Consequently, the spin correlation 
length is not minimal where the spin gap is maximal but be- 
fore. These observations support that the physical onset of 
incommensurability is the one associated with q since it is as- 
sociated with a singular behavior in ^spin, as for the MG point 
of a single chain. 



2. Spin gap and correlation lengths for fixed J2, increasing J± 

We now turn to the behavior of the spin gap when cross- 
ing the transition line at a fixed J2 = 0.6 Ji by increasing the 
transverse coupling J±. The numerical results are displayed 
on Fig. 8. Although the spin gap always remains finite, it first 
decreases before increasing at large J±. These behaviors are 
easily understood in two limiting cases : (i) at small J±, the 
effect of the transverse coupling is essentially to destabilize 
the order parameter of the dimerized chains by RVB fluctua- 
tions and concomitantly slightly increases the spin gap, (ii) at 
large J±, the spin gap is essentially controlled by the creation 
of a magnon on a rung, which increases with J± as found in 
the strong-coupling expression (14). In between, there is actu- 
ally a non-trivial local minimum found close to the transition 
but not exactly at it, and which is realized for a slightly larger 
value of J± . 

As regards the correlation lengths in the system, we exhib- 
ited previously two characteristic length scales ^spin and ^dimei- 
associated with the spin and dimer fluctuations respectively. 
We do expect that these two length-scales behave very differ- 




FIG. 8. (color online). Evolution of the spin gap and typical length 
scales with J± and a fixed J2 — 0.6 Ji. The vertical line material- 
izes the critical point. 




FIG. 9. (color online). Evolution of the triplet gap with J2 and J± 
{L = 64 sites). 



and the critical line J2 < J2.C1 J± = have a zero gap. 



ently along this cut since we know that ^dimei must diverge at 
the transition while ^spin, which is related to the inverse spin 
gap, must remain finite across the transition. This strong dif- 
ference is illustrated on Fig. 8 from which we observe nearly 
an order of magnitude between the two lengths around the 
critical point. Notice that £,^1^, as the spin gap, displays a 
minimum close to the transition but for a value of J± even 
larger than the one corresponding to the minimum of the spin 
gap- 



3. Evolution of the spin gap throughout the phase diagram 

From the two preceding cuts in the phase diagram, we ob- 
serve that the interplay between transverse coupling and frus- 
tration results in non-monotonic behaviors of the spin gap 
throughout the phase diagram. How the local minimum ob- 
served close to the transition line interpolates with the local 
maximum arising in the vicinity of the onset of incommen- 
surability in the dispersion relation is a challenging question. 
We give in Fig. 9 the typical evolution of the spin gap as a 
function of both J2 and J± for a fixed size L = 64. The line 
of local maxima survives down to the weak-coupling regime 
where it meets the line of minima at a saddle point which is 
qualitatively close to the crossing point between the transition 
and incommensurability lines of the phase diagram of Fig. 2. 
Another feature that is qualitatively well reproduced is the lin- 
ear opening of the gap with J± (forgetting logarithmic correc- 
tions recalled in Sec. I) starting from the critical phase below 
J2,c- This region of very small gaps is demanding numerically 
and the exact behavior of the gap there is beyond the scope of 
this paper Finally, the tendency at large J2 is a decrease of 
the gap for all J±. Strictly speaking, in the limit Ji = 0, 
as one recovers the situation of two decoupled ladders which 
couplings are J2 and J±, the gap is finite and starts to grow 
linearly with J±/J2- Only the corner point at J^ = Ji = 



C. Onset of incommensurability 

The onset and signatures of incommensurability is often a 
puzzling issue in frustrated models. As already introduced 
in Eq. (4), we call q the wave-vector of the real-space spin- 
correlators and q* the wave-vector at which the dispersion re- 
lation has its minimum. Having in mind that the spin structure 
factor S{k) is accessible experimentally, we can further intro- 
duce the wave-vector q at which S{k) is maximum. In unfrus- 
trated antiferromagnets, the tendency to Neel ordering usually 
translates in the fact that q = q* = q = tt. In the presence of 
frustration, we expect that, in general, the three wave-vectors 
can be different and that the value of the couplings at which 
they start to deviate from tt are not the same. For instance, 
in a gapped system, the difference between q and q is easily 
understood from a Rayleigh criterion. Indeed, by taking the 
Fourier transform of short-range correlations at wave-vector 
q, one gets a double-peak structure for S{k), which width is 
governed by S^^J^. We thus expect q to be larger than q, quali- 
tatively q = g+const. x S^~J^. Clearly, the onsetof q will occur 
after that of q (when increasing frustration). Few studies have 
investigated the comparison between q and q* among which 
we find Ref. 25. 

Typical cuts at fixed J± and increasing J2 are gathered in 
Fig. 10, showing that the larger J±, the larger the discrepancy 
between the onset and the behavior of q and q. As discussed 
previously, all wave-vectors reach 7r/2 in the large- J2 Hmit. 
The black line recalls the prediction for q* from Eq. (11) for 
comparison. In particular, we observe that q is very close to 
q* in the strong-coupling regime. 

These observations can be explained by analytical calcula- 
tions in the strong-coupling limit using a scenario proposed 
by Nomura-*^ for the onset of incommensurability. Since the 
one-magnon dispersion relation is in the ky = tt sector, we 
only consider the structure factor in this sector which ampli- 
tude is much bigger than the one at ky — 0. We define the 



o 
o 

> 

I 



3n/4 



7l/2 



1 — p — 1 1 1 1 i[*l r-< 


*-i — 1 — 1 1 1 


- q*=acos(-Jj/4J2) 


1 


A 


q (J^^i) 


\ \ I* 




A q(J^=3) 


\ \ \ 


L 


■ q(Ji=5) 


\ V \ \' 


• q(J^=io) 


\ \ \ \ 


A q(J^=3) 


- \\\\ 


□ q(J^=5) 
o q(Jj^=10) 




x[V\ 


^ 


^^^\ ■ 










* — B 






,1 7^ 



0.01 



0.1 



10 



J, /J, 



FIG. 10. (color online). Incommensurate wave-vectors q and q with 
increasing frustration J2 for different J± from DMRG results (dots) 
and compared with BOMF formula (17) (straight lines). The black 
line gives the strong-coupling prediction for q* from Eq. (11). 



spin structure-factor as 

L 

S{k) = S{k,n) =Y,e^'^^{{Sl, ~ Sl,){Sl, - SI,)} . 

x=l 

Within the BOMF approximation,^^^ we obtain the following 
expression for S{k) : 

S{k) = s^[l + di cos(fc) + rf2 cos(2A:)]-5 oc l/a;(fc) (15) 

where di and d, are the two small positive parameters intro- 
duced in Sec. Ill A 2. Remarkably, since S{k) ex l/a;(fc), in- 
commensurability in S{k) occurs on the same line J, ~ Ji/4 
as for the dispersion relation, in other words q ~ q*. How- 
ever, this is only true where the BOMF approximation is valid, 
that is, in the strong inter-chain coupling regime, and small 
deviations are found for smaller J±/Ji : for J± = Ji, we 
have checked that the position of the structure factor maxi- 
mum does not correspond to q* anymore (from Fig. 5). 

As shown on Fig. 10, incommensurability in q occurs in the 
correlation function for smaller values of J2 than in the struc- 
ture factor. Indeed, oscillations in the real space should be 
connected to singularities of S{k) in the complex plane, ^^'^^ 
the maxima on the real axis being just a signature of the sin- 
gularities. As described in Ref. 39, the onset of incommen- 
surability in the correlation function can be interpreted as a 
fusion of branch points in the complex plane. In order to dis- 
cuss these singularities, the structure factor can be rewritten 
as 



Sik) 



with P{X) = 2d2X^ + diX + l-d2 



^P(cosfc) 

We set k = a + ib with a E [0, 27r] and 6 £ M, which gives 

P(cos k) =1 + di cos(a) cosh(6) + d, cos(2a) cosh(25) 
— i sin(a) sinh(6) [M, cos(a) cosh(6) + di] . 
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FIG. II. (color online). Evolution of singularities (dots) and branch 
cuts (green lines) of the BOMF structure factor with increasing frus- 
tration. Dotted lines are the lines where 9[P(cos fc)] = 0. 



On Fig. 11, we plotted the branch cuts and singularities of 
5(fc) in the complex plane (a, fo). The discriminant of P 
changes sign for d, = d^.c ; 



d2,c = 




(16) 



For d, < d2.c, P has two real roots both inferior to —1. As 
a result, S{k) has four branch points on the axis a = tt. For 
d2 = d2,c, P can be factorized and S{k) has not any branch 
cuts but only two essential singularities on the axis a = tt. For 
^2 > cf2,c, the roots of P have a non-zero imaginary part, so 
the branch points of S{k) are not on the axis a = tt. Their 
real parts are given by the solutions of : 



2d2 cos" (a) - ( 1 + rf2 ) cos^ (a) + 



dl 



= 0. 



that is a = g or 27r — g, with, after taking the strong-coupling 
approximation for di 2, the result : 



•h 



2^1^! 



(17) 



This expression of the incommensurate wave-vector of the 
real-space correlations is in very good agreement with the 
DMRG calculations of Fig. 10. In particular, we see that the 
departure from it for both q and q* = q has the same scaling 
exponent, with (7 ~ 7r—^(J2 — J2,c)^, in agreement with the 
numerical findings of Ref. 38. We also obtain the approximate 
value for the onset of incommensurability in q as 



J2 = Ji/4J± 



(18) 



showing that it vanishes at large J± in the phase diagram of 
Fig 2, while the lines corresponding to q* and q would saturate 
at J2/J1 = 1/4. To complete the discussion, we stress that the 
incommensurability emerges in S{k) only when ^2 reaches 
di/4 (ie. J2 = t/i/4). This gives a graphical interpretation 
in the complex plane for the onset of incommensurability in 
S{k) : the position of the maximum q is then given by the 
crossing of the line cos(a) cosh(&) + di/(4d2) = (appearing 
in the imaginary part of P) with the real axis. 

Looking back at Fig. 10, the situation for J± = Ji for 
which we find that q is actually close to the strong-coupling 
prediction for q* is rather misleading and counter-intuitive. 
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In fact, this qualitative agreement results from the fact that for 
this particular choice of J^, the onset of q is not very different 
from the approximate result of Eq. (18) which itself numeri- 
cally coincides with the onset of q* . The difference between 
q and q* is only seen for sufficiently large J2, close to the 
asymptotic limit 7r/2. One may wonder whether q gets locked 
to 7r/2 before the Ji = limit is attained, in a similar way as 
q remains locked to tt below the incommensurability thresh- 
old. It is however hard to discriminate numerically between 
locking at a finite J2/J1 or an asymptotic approach to 7r/2 
despite the seeming attraction by the 7r/2 line. Still, if there is 
a locking to 7r/2, we expect it to occur in the non-perturbative 
regime where J2 ^ J± ^ Ji- 



IV. THERMODYNAMICS 

In this section, we briefly address some issues related to the 
thermodynamics of the model. The magnetic susceptibility 
and specific heat vs. temperature can be computed from full 
diagonalization of small ladders, up to L = 10. In the region 
of the phase diagram of Fig. 2 where energy scales are small 
and correlation lengths very large, we expect strong finite- 
size effects. We thus prefer to discuss the large- Jx regime in 
which all gaps are sufficiently large. Still, before embarking to 
the strong-coupling regime, we can briefly comment on some 
interesting features of the thermodynamics of the columnar 
dimer phase. As we have seen from Fig. 5, there will be four 
low-lying singlet states which contributes to the specific heat 
but not to the susceptibility on finite clusters. More interest- 
ingly, right at the critical point, as a singlet gapless mode de- 
velops, we propose that the specific heat will have a linear 
behavior at low temperature (according to the Ising universal 
exponents), while the magnetic susceptibility will be expo- 
nentially suppressed with an energy scale given by the spin 
gap. The system size we have access to are two small to un- 
ambiguously demonstrate this appealing scenario. 

We now discuss the behavior of the magnetic susceptibility 
in the rung-singlet phase. As we have reasonable analytical 
expressions for the magnon dispersion relation oj{k) in the 
strong interchain-coupling regime, we can derive the suscep- 
tibility assuming independent magnons, and compare it with 
numerical results. Using the approximate statistics of hard- 
core bosons proposed in Ref. 40, one finds the following ex- 
pression for the magnetic susceptibility : 
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FIG. 12. (color online). Magnetic susceptibility x (divided by 
{gfis)^ with g the electronic dimensionless magnetic moment and 
fiB the Bohr magneton) vs. temperature T for J± — 3Ji and 
J2 — Ji/2 from different methods. The inset displays the com- 
parison at low temperatures on a log-log scale. 



the opposite high-temperature limit, we recover the Curie 
law of decoupled rungs x{/3) — /?/4 for all calculations. 
However, corrections from this Curie-law as observed in the 
full diagonalization results (finite-size effects are negligible 
at high temperature) are not quantitatively captured by the 
one-magnon dispersion relation, signaling the relevance of 
many-magnon effects mentioned above. At intermediate 
temperatures T « Ji/ks where the maximum of the suscep- 
tibility is found, the comparison is the worst, illustrating the 
fact that the value and the position of the maximum of the 
susceptibility is a non-trivial issue. "*' 

One can however gain some qualitative information about 
the effect of frustration in the strong-coupling limit. The ex- 
pression (19) can be expanded in Ji^2/J± starting from the 
decoupled rung limit: 



x(/3) 



Xo(/3)x 
1- 



1 



1 + 3e-P-'^ 



^yf(/3Jeff)-^(/3Jeff)^ 



(20) 



with xo(/3) = /3e"'^-'^/(l + ae"*^'^^) the susceptibility 
of decoupled rungs, and the effective leg coupling Jeff = 



with the one-magnon partition function z(/3) = ^/Ji + J^- The expansion reproduces well the independent 
-ij ^^, e"'^"^*'' where uj{k) is taken from Eq. (10) and magnon result down to temperatures below the maximum. 

When the temperature is too small (/? Jx < 1), the expansion 
becomes questionable and the formula fails to reproduce the 
susceptibility. The interest of this formula is that despite its 
simplicity and strong approximations, it accounts for a non- 
trivial feature of the intermediate temperature regime which 
is the reduction of the maximum with Jeff. Since Jeff is en- 
hanced by frustration, we thus expect the maximum of the 
susceptibility to be sensitive to frustration. 



the inverse temperature /3 ~ l/ksT. In Fig. 12, this ansatz 
is compared to the full diagonalization prediction and to the 
limit of decoupled rungs Ji = J2 = for a rather large 
J± = 3Ji and frustration J2 = 0.5Ji. Since the gap from 
the strong-coupling dispersion relation is exact to second 
order in J1.2/ J±, we predict the correct behavior at low 
temperature in which many-magnon effects, comprising 
both interaction and statistical effects, are negligible. On 
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FIG. 13. (color online), (a) Evolution of the susceptibility with frus- 
tration J2 for an isotropic ladder (J^ = Ji). (b) Maximum of the 
susceptibility with frustration, (c) Temperature of the maximum. 



In order to further investigate the effect of frustration in the 
RS phase, we plot in Fig. 13(a) the results of full diagonal- 
ization calculations for a fixed J± and increasing J2. Two 
physical quantities govern the evolution of the curves : at low- 
temperature the spin gap sets an energy scale and we know 
from Fig. 7 that this gap has a non-monotonous behavior with 
J2. Extracting the gap from low-temperature measurements 
can thus provide two possible values for J2 and it is conse- 
quently not sufficient to use only the low-temperature regime. 
On the scale of Fig. 13(a), the change in the susceptibility due 
to the modification of the gap is visible. The intermediate and 
high-temperature regimes are rather controlled by Jeff/J± and 
show a systematic reduction of the susceptibility. As these 
two physical quantities evolves differently with J2, the inter- 
mediate regime where the high and low temperature behaviors 
meet might have a non-trivial behavior Indeed, the position 
of the maximum of the susceptibility displays a slightly non- 
monotonic evolution with J2, qualitatively due to the evolu- 
tion of the spin-gap (finite-size effects are small close to the 
maximum). High-temperature expansions'*^ can certainly cap- 
ture the evolution of the maximum with J2 but are beyond the 
scope of this article. 



V. CONCLUSION AND DISCUSSION 

A. Summary of main results 

In this paper, the phase diagram and most physically impor- 
tant properties of the J1-J2 two-leg ladder were computed. As 
often encountered in frustrated models, non-trivial physics is 
found : reentrance of the columnar dimer phase, low-energy 
RVB fluctuations, non-monotonic evolution of the spin gap, 
two characteristic lengths, and the onset of incommensura- 
bility. Apart from its experimental relevance which will be 
discussed in the next section, this model provides a simple 



physical picture and possibly experimentally relevant exam- 
ple for the melting of a dimer crystal. Furthermore, the sur- 
prising properties of the rung-singlet phase in the presence of 
frustration can be understood analytically and quantitatively 
in the large- J^ regime. In particular, it provides a clear exam- 
ple and quantitative predictions on the difference between the 
different wave-vectors (denoted q, q and q* in the manuscript) 
signaling incommensurability. We eventually argue that the 
temperature-dependence of the magnetic susceptibility can be 
well understood in the low and large-T limits, but that the 
intermediate regime where the maximum is found and its evo- 
lution with frustration cannot be quantitatively understood by 
approximate methods. 



B. Experimental relevance to BiCu2P06 

Frustration in ID and quasi- ID quantum magnets plays a 
fundamental role. A famous example is the spin-Peierls tran- 
sition observed almost two decades ago in the inorganic mate- 
rial CuGeOs,'*^^ believed to be a good realization of a Ji — J2 
chain with a frustrated ratio J2/J1 ~ 0.35.** More recently, 
a new class of frustrated ferromagnetic spin chains which 
display spiral (incommensurate) correlations and ferroelectric 
behavior, including LiCuV04, LiCu202 and Li2ZrCu04 has 
been intensively studied.'*^ 

Here, we want to focus on the spin ladder system 
BiCu2P06^''^'*^ which is believed to be a realization of the 
frustrated ladder model studied in this paper. So far, most of 
the studies agree on the fact that BiCu2P06 has a ladder-like 
structure with a second-neighbor antifeiTomagnetic exchange 
along the leg direction. '^"^'''*^ However, various studies dis- 
agree regarding the amount of frustration J2/J1, the strength 
of the rung exchange J±, and the inter-ladder couplings. In- 
deed, the size of the spin gap Ag/Ji ^ 0.3 cannot be sim- 
ply explained by an isolated frustrated ladder model where it 
is overestimated.^' Additional 2D (and 3D) couplings'*^-'*^ be- 
tween ladders are expected to play an important role: (i) to 
reduce the spin gap and (ii) to account for impurity-induced 
ordering at finite temperature.'*'^ Another argument in favor 
of non-negligible inter-chain effects comes from the field- 
induced magnetization curve. ^' Taking all these facts into ac- 
count, it is reasonable to argue that the most realistic model for 
the spin ladder material BiCu2P06 should integrate both in- 
gredients: in-leg frustration J2 and higher-dimensional inter- 
ladder couplings. However, depending on the properties we 
want to investigate, dimensionality effect may matter or not. 
For instance, impurity-induced ordering at finite temperature 
is clearly a 3D effect,'*^ as well as NMR lines broadening at 
low temperature."*^ On the contrary, the incommensurate re- 
sponse is expected to be a ID effect simply because the frus- 
tration is present only along the legs. Looking at Fig. 10, 
one can see that if Ji ~ J± (which seems to be the case for 
BiCu2P06), an incommensurate response is already expected 
for J2/J1 ^ 0.25. Inelastic Neutron Scattering experiments 
on single crystal samples may allow for a direct evaluation of 
the incommensurate wave-vectors q and q* which should help 
to determine the frustrated character of this material. 
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The evaluation of the amount of frustration in BiCu2P06 
from susceptibility measurement requires to fit the full range 
of temperature to capture both the spin gap and the max- 
imum Xmax of the susceptibility on equal footing. Inter- 
ladder interactions will also affect the spin gap and the many- 
magnons effects relevant in the intermediate and high temper- 
ature regimes. A simple unfrustrated model of coupled lad- 
ders has been shown using QMC simulations (unable to tackle 
frustrated models because of the sign problem) to correctly 
reproduce the low temperature part of the susceptibility while 
the maximum of x was not correctly described. "^^ Conversely, 
the high temperature part and the reduced value of Xmax was 
correctly captured by a ID model including frustration (using 



ED) but at the expense of a spin gap largely overestimated.^' 
Therefore, we believe that a precise experimental determina- 
tion of the incommensurate response will be of crucial interest 
in order to estimate the amount of ID frustration in this sys- 
tem. 
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